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ABSTRACT 

We calculate NLTE line-driven wind models of selected O stars in the spectral range of 04 
to 09 in the Small Magellanic Cloud (SMC). We compare predicted basic wind properties, 
i.e. the terminal velocity and the mass-loss rate with values derived from observation. We 
found relatively good agreement between theoretical and observed terminal velocities. On the 
other hand, predicted mass-loss rates and mass-loss rates derived from observation are in a 
good agreement only for higher mass-loss rates. Theoretical mass-loss rates lower than ap- 
proximately 10~ 7 M©year -1 are significantly higher than those derived from observation. 
These results confirm the previously reported problem of weak winds, since our calculated 
mass-loss rates are in a fair agreement with predictions of Vink et al. (2001). We study multi- 
component models for these winds. For this purpose we develop a more detailed description 
of wind decoupling. We show that the instability connected with the decoupling of individual 
wind elements may occur for low-density winds. In the case of winds with very low observed 
mass-loss rate the multicomponent effects are important for the wind structure, however this 
is not able to consistently explain the difference between predicted mass-loss rate and mass- 
loss rate derived from observation for these stars. Similar to previous studies, we found the 
dependence of wind parameters on the metallicity. We conclude that the wind mass-loss rate 
significantly increases with metallicity as M ~ Z - 67 , whereas the wind terminal velocity on 
average depends on metallicity only slightly, namely Voo ~ z o oe (for studied stars). 

Key words: stars: winds, outflows - stars: mass-loss - stars: early-type - hydrodynamics - 
instabilities - galaxies: Magellanic Clouds 



1 INTRODUCTION 

In recent few years 8m class telescopes became routinely available 
for the stellar research. This enabled detailed study of many stars 
in the Local Group. Clearly, many astrophysically important types 
of stars can be studied from another perspective. This is especially 
true for stars from the Small Magellanic Cloud (SMC). Their lower 
metallicity compared to the solar value (the average metallicity of 
in dividu al elements is Z/Zq — 0.2, e.g. Venn 1999, Bouret at 
al. l2003l) enables to study the stellar properties and evolution with 
respect to metallicity. 

One of the most important properties of hot stars, that can in- 
fluence both the observed spectrum and the stellar evolution, is the 
presence of the stellar wind. The existence of significant depen- 
dence of the basic wind properties (i.e. the amount of mass expelled 
from the star per unit of time, the mass-loss rate and the wind ve- 
locity at large distances from the star, the terminal velocity) on the 
metallicity has been anticipated already at the very beginning of the 
theoretical study of hot star winds (Abbott 1982, Kudritzki et al. 
Il987t) . Although the most abundant elements in the observed uni- 
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verse, hydrogen and helium mostly contribute to the wind density, 
their contribution to the radiative acceleration is very small (e.g. 
Abbott fl982h . Heavier elements like carbon, nitrogen, oxygen or 
iron are much more important for the wind acceleration of present 
hot stars. This is because heavier elements have effectively much 
more lines that are available to absorb the stellar radiation and to 
accelerate the stellar wind. Clearly, in most cases the radiative force 
shall be higher for higher metallicity. Consequently, the basic wind 
properties shall depend on the metallicity. 

Although relatively simple expressions for the dependence of 
the mass-loss rate and the terminal velocity on the metallicity can 
be obtained even on the basis of the standard CAK (Castor, Abbott 
& Klein 1 19751) theory (e.g. Puis et al. ll998h . th e deta iled depen- 
dence is probably more complicated (Puis et al. 1998, Vink et al. 
2001). It is clear that more precise (NLTE) wind models are neces- 
sary to study this problem in detail. 

There are several independent NLTE wind models that can 
be used for the study of the metallicity dependence of the wind 
properties. These codes differ in the level of sophistication of the 
treatment of the wind problem. One of the most important aspects 
that influences the reliability of individual NLTE wind models is 
the inclusion of proper bound-free and bound-bound blanketing in 
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the UV domain. While models of Vink et al. 1 200 1 , hereafter VKL) 
use Monte Carlo method for the solution of the radiative transfer 
equation, models of Pauldrach et al. ( 2001) s olve t his equations in 
detail. Wind models of Grafener & Hamann (2005), that were suc- 
cessfully used for the calculation of wind models of WR stars, use 
comoving-frame formulation of the radiative transfer equation. All 
these models to some extent properly take into account blanketing 
effects in the UV domain that are important for the correct calcu- 
lation of the ionization balance and of the radiative force. Hence, 
these models are able to reliably predict the most important wind 
parameter - the mass-loss rate and its metallicity dependence. 

However, also the wind velocity field and especially the ter- 
minal velocity can be influenced by the metallicity. To study this 
effect, the solution of hydrodynamic equations is necessary. Since 
models of VKL, that are widely used in hot star evolutionary cal- 
culations, do not solve these equations, these models assume a pre- 
specified velocity law and cannot be used for the predictions of the 
velocity structure. This is an important difference between models 
of VK L and e.g. Pauldrach et al. 1200 ll) or Grafener & Hamann 
(2005). Hydrodynamical win d mod els for different metallicities 
were calculated by Kudritzki (2002). According to these models 
the wind terminal velocity decreases with decreasing metallicity, in 
agreement with SMC observations (see an extensive compilation of 
Kudritzki & Puls l200d) . 

However, the wind metallicity does not vary the radiative force 
only. The momentum acquired by the heavier elements from the ra- 
diative field is transferred to the bulk wind component (i.e. hydro- 
gen and helium) via the Coulomb collisions. Since the frictional 
force caused by these collisions depends on number densities of 
heavier elements and hydrogen-helium component, the effects in- 
duced by the collisions are more important in the low-metallicity 
environment (Kudritzki l2002l Krticka et al. 120031) . NLTE models 
available to study this problem were presented by Krticka & Kubat 
<2004 hereafter KK1). 

Some hot stars with low luminosities (and consequently low 
mass loss rates) show signatures of very weak winds, much weaker 
than that ded uced from st andard wind theory (Bouret et al. l2003t 
Martins et al. l2004l : |2005t) . Sinc e this discrepancy is based mostly 
on wind models of Vink et al. (2001), it would be interesting to 
test this result also with independent NLTE wind models and to test 
whether this discrepancy can be caused by multicomponent effects. 

The study of the metallicity dependence of stellar wind prop- 
erties is not important only for our understanding of stellar winds 
or massive stars itself. Since stellar winds play a role in the chem- 
ical enrichment of galaxies and are responsible for input of a large 
amount of momentum and energy in the interstellar medium, the 
detailed knowledge of variations of basic wind parameters with 
metallicity is important also for other fields of astrophysics. 

As we have discussed, wind properties of hot stars shall de- 
pend on the stellar metallicity. Low-metallicity environment of 
SMC offers a unique possibility to test these available predictions. 
To do so, we present here wind models for selected SMC hot stars. 



2 MODEL DESCRIPTION 

The models applied in this paper were in detail described by KK1. 
Here we only summarise the basic model properties and refer an 
interested reader to this paper. 

Models assume spherically symmetric and stationary stellar 
wind. Occupation numbers of selected atoms and ions (see Ta- 
ble0 are obtained by the solution of statistical equilibrium (NLTE) 



Table 1. Atoms and ions included in the NLTE calculations. In this table, 
level means either an individual level or a set of levels merged into a super- 
level. 



Ion Levels 



Ion Levels 



Ion Levels 



Hi 


9 


NelV 


12 


SIV 


18 


Hii 


1 


Ne V 


17 


Sv 


14 


He I 


14 


Ne VI 


1 


S vi 


16 


Hell 


14 


Nail 


13 


S vii 


1 


Hem 


1 


Nam 


14 


Arm 


25 


CII 


14 


NalV 


18 


Ariv 


19 


Cm 


23 


NaV 


16 


ArV 


16 


Civ 


25 


NaVI 


1 


Ar VI 


1 


Cv 


1 


Mgll 


14 


Call 


16 


Nil 


14 


Mgm 


14 


Cam 


14 


Nm 


32 


Mgiv 


14 


Caiv 


20 


Niv 


23 


Mg V 


13 


Cav 


22 


Nv 


13 


Mg VI 


1 


Ca VI 


1 


Nvi 


1 


Aim 


14 


Fein 


29 


Oil 


50 


Aliv 


14 


Feiv 


32 


Om 


29 


Alv 


1 


Fev 


30 


OIV 


39 


Sim 


12 


Fe vi 


27 


OV 


14 


Si IV 


13 


Fe VII 


1 


O VI 


20 


Siv 


15 


Nim 


36 


O vii 


1 


Si VI 


1 


NilV 


38 


Nell 


15 


So 


14 


Niv 


48 


Nelll 


14 


S in 


10 


Ni VI 


1 



equations together with the radiative transfer equation. The radia- 
tive transfer equation in lines is solved in the Sobolev approxima- 
tion (Sobolev fl947l Castor ll974h . whereas the continuum radiative 
transfer equation is solved by Feaut rier m ethod in t he sp herical co- 
ordinates (see Mihalas & Hummer 1 19741 or Kubat fl993l) . Derived 
occupation numbers are used to calculate the radiative force (in the 
Sobolev approximation) and the radiative cooling/heating t erm (u s- 
ing the thermal balance of electrons method, Kubat et al. ( 1999)). 
This enables to solve the hydrodynamic equations. These proce- 
dures are iterated to obtain consistent model structure. Finally, wind 
mass-loss rates and terminal velocities of studied stars can be ob- 
tained from our models and compared with observations. 

Although our models involve some simplifying assumptions 
(especially the simplified treatment of the radiative transfer (split- 
ting of the radiative transfer in continuum and in the lines, neglect 
of line overlaps) compared to more advanced models of e.g. Paul- 
drach et al. ( 2001 ), VKL or Grafener & Hamann ( 2005)), they were 
able to correctly predict basic wind parameters of late O stars (see 
KK1) and of A supergiant (see Krticka & Kubat 2004). Moreover, 
our models have some advantages compared to some other mod- 
els available in the literature, for example the direct calculation of 
the radiative force without using of force multipliers or multicom- 
ponent treatment of model equations (see Krticka & Kubat l200ll 
hereafter KKII). 

Presented models were only slightly modified with respect to 
the status described by KK1. First, we included accelerated lambda 
iterations in continuum (or, more precisely, approximate Newton- 
Raphson iterations) based on Rybicki & Hummer 1 19921) paper and 
Ng acceleration (Ng[l974J) to accelerate the convergence of the sys- 
temof statistical e quilibr ium equations (see also Hillier & Miller 
( 1998), or Hubeny 12003T) for a review). This will be described in 
a following paper (Krticka & Kubat, in preparation). Second, our 
previous set o f included atomic models based on TLUSTY files 
(Hubenv ll98i Hubeny & Lanz ll992L Hubeny & Lanz ll995t Lanz 
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& Hubeny [2003) or Opacity P roject fSeaton Il987l Lup & Prad- 
ha nll989l Sawey & Berrin gtonl 19921 Seaton et al. ll992l Butler et 
al. 1993, Nah ar & P radhan 1993) and Iron Project (Humm er et al. 
Il993l Baut istall99d Nahar & Pradh anll996t Zhang ) 19961 B autista 
& Pradhan fl99*T zhang & Pradhan [l93chen & Pradhan fl999T) 
data was slightly extended (see Table0. These mentioned changes 
do not significantly influence the derived results. To complete the 
list of atomic databases used, the oscillator strengths necessary for 
the calculation of the radiative force are extracted from the VALD 
database (Piskunov et al. ll995l Kupka et al. ll999h . 

The boundary radiative flux is taken from grid of line- 
blanketed plane-parallel model atmospheres OSTAR2002 (Lanz & 
Hubeny 2003) instead of from H-He spherically symmetric models 
of Kubat (2003). Since line-blanketed fluxes have generally lower 
flux in the UV region (where are many lines important for radia- 
tive driving), obtained mass-loss rates are slightly lower (roughly 
1.4 x) than that derived using H-He fluxes. However, this differ- 
ence is lower than the dispersion of mass-loss rates for Galactic O 
stars (see KK1) and does not significantly influence predicted wind 
terminal velocities (the difference between the terminal velocities 
is about lOOkms -1 ). 



3 CALCULATED WIND MODELS 
3.1 Parameters of studied SMC stars 

Since our main intention for future studies is to study weak winds 
of B stars, for the present study we selected only those cooler 
(i.e. with effective temperatures T c s < 42 000 K) O SMC stars, 
for which at least reliable estimate of their mass-loss rate is avail- 
able in the literature. We tried to omit those stars for which their 
wind parameters are uncertain and which are binaries. Moreover, 
we also aim to base the list on broader surveys with larger number 
of individual stars studied. Stellar parameters and wind parame- 
ters of these stars are given in Tab.[2] Stellar effective temperatures 
and radii arc taken from Puis ct al. 1 1996, hereafter P96). Bouret ct 
al. <2003l hereafter B03) and Massey et al. fe004l hereafter M04). 
Whereas parameters given by P96 were obtained by models without 
wind-blanketing, parameters given by B03 and M04 were derived 
using models with wind-blanketing. To study the thin-w ind prob- 
lem in detail we added also stars from Martins et al. 120041 hereafter 
Mr04) sample. These stars are suspected Vz stars, i.e. stars close to 
the ZAMS. They exhibit much lower wind spectral signature than 
that predicted from standard theory. On the other hand there are 
only upper limits of their observed mass-loss rates and lower limits 
of their observed terminal velocities available. 

SMC stars are in some sense more suitable for the test of the- 
oretical models than stars from our Galaxy. Since the distance to 
SMC is known with relatively high precission, the stellar radius 
and mass-loss rate may be also derived more reliably. 

For our study we adopted evolutionary stellar masses either 
derived by B03 or by us using evolutionary tracks of stars with 
initial metallicity Z/Zq = 0.2 calculated by Charbonnel et al. 
1 1993). The use of the evolutionary masses is a relatively important 
assumption that can significantly influence the results derived due 
to the discrepancy between stellar masses of hot stars derived from 
evolutionary tracks and from spectroscopy (Herrero et al. 1992). 
Since KK1 in their analysis used evolutionary masses, we also use 
evolutionary masses here, but in fact for many stars from KK1 sam- 
ple these masses are nearly equal. Our SMC sample is not homo- 
geneous in that sense, since for some stars these masses are equal, 



but for a significant part of stars the evolutionary mass is roughly 
1.5x higher than the spectroscopic one. The use of spectroscopic 
masses instead of the evolutionary ones would help to obtain a bet- 
ter agreement between observation and theory for some stars, how- 
ever it would cause differences for some others. 

When available, the abundances of individual elements were 
taken from the literature, however in other cases we assum ed av- 
erage value Z/Zq = 0.2 derived for SMC stars (e.g. Venn ll999L 
M04). We use Galactic helium abundance. 

There are indications that stellar winds of O stars are clumped 
(e.g. B03). From the observational point of view, the possible wind 
clumping decreases the wind mass-loss rate inferred from the ob- 
servation because the line profiles of clumped wind mimic those 
with higher mass-loss rate. According to the numerical simulations 
of wind instability (e.g. Feldmeier et al. ll997l Runacres & Owocki 
2002) the theoretical mean mass-loss rate is nearly the same for 
smooth and structured winds. Thus, if the observations really show 
signatures of clumping at all spectral regions where the stellar wind 
is observed, then, according to our present knowledge, the theo- 
retically predicted values of mass-loss rates should basically cor- 
respond to the values derived from observation with account of 
clumping. For our study we adopted the values derived from obser- 
vations assuming "smooth" winds partly because these values were 
for larger SMC sample (to our knowledge) derived only by B03 and 
partly because the models of wind clumping are still schematic. 

Wind parameters adopted for the comparison with theoretical 
values will be discussed individually for those stars, for which it 
is necessary. Note that in the text we will term the mass-loss rate 
estimated from observation as "observed mass-loss rate", although 
one should keep in mind that this quantity cannot be directly de- 
rived from spectra and that it is model-dependent. 

NGC 346 WB 1 T his is a multiple system (Heydari-Malayeri & 
Hutsemekers 1 1 99 ll) . The terminal velocity derived for this system 
by P96 «oo = 2650 km s -1 is marked as uncertain due to the 
complexity of the absorption profile. P96 note that another possi- 
ble value of the terminal velocity is Voc — 2250 km s^ 1 . Prinja & 
Crowther 1 1998) obtained value of edge velocity (velocity at which 
the line profile meets the continuum level) as u c d gc = 2830 km s~ 1 
for Nv lines. Using their derived approximate relation v x = 
0.8i> e dge we obtain terminal velocity Voo = 2260 km s -1 . Thus, 
we adopted = 2250 krns -1 as a value of the terminal velocity. 

NGC 346 WB 4 The observed terminal velocity by P96 v x = 
1550 km s -1 is quoted as uncertain. The typical terminal veloc- 
ity for stars of similar spectral type is higher, the terminal ve- 
locity calculated as twice the escape velocity (roughly suitable 
for SMC stars, B03) is 1900kms _1 . Using approximate relation 
Woo = 0.8w e dge of Prinja & Crowther 1 19981) and their measure- 
ment of edge velocity of C IV lines we obtain = 1950 km s _1 . 
Hence, we adopt this value of terminal velocity. 

NGC 346 WB 6 The edge velocity obtained for this star by Prinja 
& Crowther 1 19981) w e d gc = 1925kms _1 is lower than the ter- 
minal velocity Vaa = 2250 km s^ 1 derived by P96 or v x — 
2300 km s -1 by B03. This difference probably illustrates the fact 
that the correct determination of wind velocities for SMC is diffi- 
cult due to their low wind density. 

This star was independently studied by P96, who obtained 
slightly lower value of the effective temperature than B03 (T e g = 
40 000 K). 
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Table 2. Stellar and wind parameters of selected SMC stars. Stellar parameters (radius R t , the effective temperature T cf j and th e meta llicity relative to the 
solar value Z /Z^) were adopted either from Puis et al. <199d hereafter P96), Bouret et al. 120031 hereafter B03), Massey et al. j2004, hereafter M04) and 
Martins et al. 1 2004, hereafter Mr04). For the stellar mass M we assume values derived using evolutionary tracks. For stars with metallicities denoted as B03 
and Mr04 we adopted detailed abundance determinations from B03 and Mr04 for C, N, O, Si, S and Fe and Z/Zq =0.2 for other heavier elements. Observed 
wind parameters (i.e. the mass-loss rates M and the wind terminal velocities Voc) were also mostly taken from P96, B03, M04 and Mr04 (when available), 
however see the discussion for individual stars. Predicted values of wind parameters were derived by our code. For stars for which v^o value is not given in 
the table authors provide only lower limits that is in agreement with our predicted value of Voo ■ 



Star 




Stellar parameters 




Mass loss rates M 


Terminal velocities Doc 


Source 




Sp. 




M 




Z/Zq 


observed 


predicted 


observed 


predicted 








[R Q ] 


[Mq] 


[K] 




[M yr" 1 ] 


[M Q yr" 1 ] 


[kms ] 


[kms ] 




NGC 346 WB 1 


04III 


23.3 


95 


42 000 


0.2 


4.8 10~ 6 


3.5 10" 6 


2250 


2240 


P96 


NGC 346 WB 4 


05.5V 


14.2 


53 


42 000 


0.2 


< 1 10~ 7 


7.8 10~ 7 


1950 


2200 


P96 


NGC 346 WB 6 


04V 


11.2 


40 


41 500 


B03 


2.7 10~ 7 


3.5 10~ 7 


2300 


2180 


B03 


NGC 346 MPG 368 


05.5V 


10.6 


38 


40 000 


B03 


1.5 10~ 7 


1.810- 7 


2100 


2580 


B03 


NGC 346 MPG 113 


06V 


7.8 


33 


40 000 


B03 


310" 9 


4.8 10~ 8 




3370 


B03 


N81 #2 


06.5 


7.9 


31 


40 000 


Mr04 


<10-8 


410-8 




2880 


Mr04 


AzV75 


05III 


25.4 


92 


40 000 


0.2 


3.5 10~ 6 


3.8 10- 6 


2100 


2120 


M04 


N81 #1 


07 


10.3 


34 


38 500 


Mr04 


<io~ 8 


910- 8 




2530 


Mr04 


AzV 26 


061 


27.5 


86 


38 000 


0.2 


2.5 10" 6 


3.8 10~ 6 


2150 


1850 


M04 


AzV 232 


071a 


29.3 


93 


37 500 


0.2 


5.5 10~ 6 


4.1 10- 6 


1400 


1880 


P96 


AzV 207 


07V 


11.0 


33 


37 000 


0.2 


110~ 7 


1.0 10- 7 


2000 


2060 


M04 


N81 #11 


07.5 


6.9 


23 


37 000 


Mr04 


<io~ 9 


210-8 




2270 


Mr04 


N81 #3 


08.5 


5.0 


19 


36 000 


Mr04 


<310~ 9 


310- 9 




2260 


Mr04 


AzV 238 


09III 


15.5 


37 


35 000 


0.2 


T.310" 7 


2.8 10~ 7 


1200 


1830 


P96 


NGC 346 MPG 487 


06.5V 


10.2 


25 


35 000 


B03 


310~ 9 


3.010- 8 




2300 


B03 


AzV 469 


08.511 


21.2 


38 


32 000 


0.2 


1.8 10~ 6 


4.0 10- 7 


2000 


2010 


M04 


NGC 346 MPG 12 


09.7V 


10.1 


21 


31000 


B03 


HO" 10 


1.710-8 




1820 


B03 



AzV 232 Crowther et al. (2002) applied models with wind blan- 
keting to study the stellar parameters of this star and obtained 
much lower value of the effective temperature than P96 (T c g = 
32 000 K). However, we calculated wind parameters with this new 
determination of stellar parameters (i.e. the effective temperature, 
stellar mass, radius and abundances of C, N, O) and we obtained 
too low value of the mass loss rate (Af = 9 10~ 7 Mq yr- 1 ), much 
lower than the observed value. Also models of VKL predict with 
this new parameters about five times lower mass-loss rate than the 
observed value. Probably, this may be caused by the fact that this 
star has a very peculiar chemical composition (Z(C) /Zq(C) = 
0.07, Z(0)/ Zp,(0) = 0.1, whereas Z(N)/Z Q (N) = 2.0, 
Crowther et al. 2002) and abundances of other elements which are 
also important for driving of the stellar wind were not determined. 
Thus, to keep our sample more homogeneous, we used the stellar 
parameters derived by P96. 

3.2 Comparison of calculated wind parameters with 
observation 

In Fig.Q we compare calculated wind terminal velocities for se- 
lected stars with observed values. While for some stars there is 
a very good agreement between observed and calculated terminal 
velocities, for some stars the agreement is worse. Before we will 
discuss this problem in a next section let us conclude, that derived 
scatter between observed and calculated wind terminal velocities is 
slightly higher than that for Galactic stars (see KK1). 

Generally, it is difficult to assess the accuracy of determina- 
tion of wind parameters. P96 state that due to the edge variability 
and contamination of the edges of wind absorption lines by under- 
lying photospheric lines the terminal velocities should be regarded 
as accurate up to ±10 %. Taking this as a typical error of terminal 
velocity determination we conclude, that with exception of stars 
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Figure 1. Comparison of calculated and observed terminal velocities. Stars 
for which there is no reliable estimate of their terminal velocities (in fact 
those stars that exhibit very weak winds) were excluded from the plot. Line 
denotes one to one relation. 



NGC 346 MPG 368, AzV 26, AzV 232 and AzV 238 the rest of 
terminal velocity estimations lie within the mentioned uncertainty 
interval. 

According to the predictions of stellar wind theory the wind 
terminal velocity is correlated with the surface escape veloc- 
ity v esc . For Galactic O stars the value of Woo/ifec does not depend 
on their effect ive tem perature and has a mean value of about 2.5 
(Lamers et al. Il995t) . Thus, we also plot the ratio of v^/v^ for 
studied stars (see Fig. [2] where we compare ratio of Vac/v^ de- 
rived using theoretical and observed values of t>oo). The mean value 
of i>oo/«esc ~ 2.3 is slightly lower than that obtained by KK1 for 
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Figure 2. Comparison of the ratio of the wind terminal velocity Vrx, to the 
surface escape velocity u esc derived using theoretical and observed values 
of Vac ■ Note that the observed values for stars with only lower limit of the 
terminal velocity available (missing Vac in Tab. l2l are not plotted in this 
graph. 
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Figure 3. Comparison of calculated mass-loss rates (either by us or by 
VKL) and mass-loss rates derived from observation. Points corresponding 
to the upper limits were obtained using our model predictions and upper 
limits derived from observations. Line denotes one to one relation. 



Galactic stars. Our calculated Voo /v esc ratio is slightly lower than 
recent observational finding of Evans et al. 1 2004), who derived the 
median Uoo/fesc = 2.63 for SMC stars with effective temperatures 
higher than 24 000 K. Since B03 obtained for their (however lim- 
ited) sample Vac/v^ ~ 2.3 we conclude that our calculations are 
in agreement with other studies and that our results may indicate 
that SMC hot stars terminal velocities are slightly lower than that 
of Galactic hot stars. 

Similarly to Galactic O stars, there is a large scatter of both 
observed and theoretical v x /v ex values (see Fig. |2j- The scatter 
of theoretical values of Woo/^esc is slightly higher for SMC stars 
than for Galactic stars (KK1). The origin of this scatter is prob- 
ably partly connected with high sensitivity of terminal velocities 
on detailed wind parameters in the outer wind regions (see Puis et 



al. l200d and also Sect. 14. 3t and also with uncertain values of the 
terminal velocities. 

The comparison of calculated mass-loss rates and mass- 
loss rates derived from observations in Fig. [5] shows relatively 
good agreement for stars with higher mass-loss rates (M > 
10~ 7 Mq yr -1 ), in many cases better than for stars in our Galaxy 
(KK1). This is probably due to the fact that the distance to SMC 
is known with a relatively high precision and, thus, the basic stel- 
lar parameters are known also with a high precision (probably with 
exception of some systematic effects like wind blanketing effect 
which may influence parameters of both stellar groups). However, 
there is a significant disagreement between theoretical and ob- 
served values for lower mass-loss rates (M < 10~ 7 M0 yr _1 ). 
In this case the predicted mass-loss rates are more than ten times 
higher than that derived from spectral analysis of observed data. 
This is not only a problem of our models, since also predictions 
of VK L show the same behaviour (see B03; Martins et al. l2004t 
2005). To demonstrate this conclusion, we have added the mass- 
loss rate predictions of studied stars calculated using VKL recipe 
into Fig.[3] The possible origin of this discrepancy will be discussed 
in Sect.|5| 

Our results for possible young Vz stars with thin winds (M < 
1O _7 M0 yr -1 ) from Mr04 sample do not show that the system- 
atic disagreement between predicted mass-loss rates and mass-loss 
rates derived from observation is higher than for generally older 
stars with thin winds from B03 sample (see Fig.|5J. 

For stars with very small mass-loss rates (M < 
lO _7 M0yr _1 ) only upper limits of their terminal velocities are 
available in the literature. Predicted terminal velocities of these 
stars are consistent with these upper limits. Note however that this 
is not a check of the models because the standard wind theory pre- 
dicts that the terminal velocities and mass-loss rates are related - if 
the mass-loss rate of these stars is really much lower than the pre- 
diction derived in this paper, then their terminal velocities may be 
much higher. 

The mass-loss rates derived by VKL recipe are slightly higher 
than that of us although KK1 found relatively good agreement be- 
tween these theoretical rates. This is caused by a different boundary 
flux used. The line-blanketed fluxes from OSTAR2002 grid used in 
this work have lower flux in the UV domain (this domain is im- 
portant for the line-driving) than H-He models used by KK1, con- 
sequently derived mass-loss rates are slightly lower. However, the 
effect of different boundary fluxes on terminal velocities is small. 



3.3 Wind momentum luminosity relationship 

Wind momentum-luminosity relationship (see Kudritzki & Puis 
fcoooh and references therein) may provide an independent method 
for the determination of stellar and consequently also galactic dis- 
tances. However, to achieve this, detailed calibration of this re- 
lationship is necessary, especially with respect to the metallicity. 
Low-metallicity environment of SMC provides an ideal tool for 
such a calibration. 

We compare theoretical modified wind momentum-luminosity 
relationship M«oo (-R»/Rq) obtained for considered stars us- 
ing our NLTE wind models with a relationship derived from ob- 
servations of these stars (see Fig. |4j. We excluded observed val- 
ues for stars for which only upper limits of their terminal veloc- 
ities were derived from observation. We conclude that there is a 
relatively good agreement between calculated values and those de- 
rived from observation. Note however that the agreement for ex- 
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Figure 4. Comparison of calculated (crosses) and observed (circles) mod- 
ified wind momentum for considered stars. Note that observed values for 
stars with unknown terminal velocity were excluded from this plot since it 
is possible to obtain at most only lower limits to their modified wind mo- 
mentum (for stars with known observed mass-loss rates). Linear regression 
of calculated data of SMC and Galactic (KK1) stars is also plotted. 



Table 3. Comparison of modified wind momentum-luminosity relationship 
(see Eq. Q) taken from different sources. SMC denotes sample considered 
in this paper, either calculated from theoretical values or from values de- 
rived from observation (again excluding the observed values for stars with 
very weak winds). 



Sample 


log D 


X 


Galactic (Vink et al. 20001 


18.68 ±0.26 


1.826 ±0.044 


Galactic (KK1) 


18.7 ±2.3 


1.83 ±0.40 


SMC (theoretical) 


16.6 ±0.2 


2.08 ± 0.04 


SMC (observed) 


16.6 ± 2 


2.1 ±0.4 



eluded stars with unknown terminal velocity (for stars with M < 
10 -7 Mq year -1 ) is very poor (using either lower limits of their 
terminal velocities or terminal velocities derived from our calcula- 
tions), since their predicted mass-loss rates are much higher than 
the observed ones. On the other hand standard theory predicts that 
terminal velocities and mass-loss rates are related. In such a case 
the terminal velocities of these stars could be much higher and 
the agreement between predicted and observed wind momentum- 
luminosity relationship could be improved. 

We have also calculated the linear regression of both theoret- 
ical and observed modified wind momentum-luminosity relation- 
ship for considered stars 



log Mvoo (R,/K 



,1/2 



Q j ■ j =xlog(L/L )+lo gj D o , (COS) 

(1) 

and compared it with the theoretical values derived for Galactic 
stars (Tab.|3J- First, due to relatively good agreement of mass-loss 
rates there is a good agreement between theoretical and observed 
modified wind momentum-luminosity relationship for SMC stars 
(exluding stars with thin wind). Moreover, due to lower mass-loss 
rates of SMC stars the Do value is significantly lower than that 
for Galactic stars. Also the slope is slightly different. Roughly 5 — 
6x lower modified wind momentum derived in the present study 
compared to that derived by KK1 for Galactic stars is mostly due 
to lower wind mass-loss rates caused by low metallicity of SMC 



and to a lesser extent by a downward revision of mass-loss rates 
due to the use of blanketed model atmospheres. 



4 INFLUENCE OF THE UNCERTAINTIES OF STELLAR 
PARAMETERS DETERMINATION 

Individual stellar parameters (like the effective temperature, mass 
or radius) are not in some cases determined with a high degree of 
precision. This is especially true for hot stars for which large uncer- 
tainties of derived stellar radii (for Galactic stars due to uncertain- 
ties of the determination of their distance) and stellar masses (see 
Herrero et al. ll99l Lanz et al. ll996l) may exist. 

The uncertainties of determination of the stellar parameters in- 
fluence also the predicted wind parameters (i.e. the mass-loss rates 
and the terminal velocities). To describe this in detail, we perform 
study of variations of predicted wind parameters with small change 
of stellar para meters. From the scaling of modified CAK theory 
(Kudritzki et al. 1989) 



M~ N 1/a 'L 1/a ' [M(l-r)] 1 



-l/o 



(2) 



where iVo is connected with number of lines that effectively drive 
the stellar wind, F is the Eddington parameter and 

a' — a — 5, (3) 

where a and 5 are usual CAK parameters, so a' is typically equal 
to about 0.5. Hot star wind mass-loss rate mostly depends on the 
stellar luminosity L. With decreasing metallicity parameter A^o 
decreases causing a decrease of mass loss rate. Finally, with de- 
creasing effective mass M (1 — F) mass loss rate increases. Since 
Eq. does not take into account the variations of a' with e.g. the 
effective temperature, we also include the scaling of Vink et al. 
(2000) obtained for Galactic OB stars 



M 



3 T e ff («oo/fesc)~ 



with result of metallicity dependence derived by VKL as 

M ~ Z 69 . 



(4) 



(5) 



From the classical CAK theory it is known that the terminal veloc- 
ity V pg depe nds mostly on the escape velocity v esc (see Kudritzki & 
Puis |2000) for a review). 

Here we do not aim to fit expressions like Eq. J2j for large 
amount of stars, however just to study variations of wind parame- 
ters with varying stellar parameters. 



4.1 Effective temperature and stellar mass 

Comparison of mass-loss rates and terminal velocities for studied 
stars calculated with original set of parameters and with stellar ef- 
fective temperature by 1000 K higher is given in Fig. [5] Appar- 
ently, for higher effective temperature the mass-loss rate is higher, 
in agreement with other theoretical predictions (see Eqs. J4j). 
From our calculations it is possible to derive that the mass-loss rate 
scales with the effective temperature as 



M 



rrilO.64 

Jeff : 



(6) 



from which the effective parameter a 1 — 0.38. This sc aling is in 
a relatively good agreement with scaling of Vink et al. 120001 see 
Eq. J4j), which implies M ~ Tjf 8 . The dependence of the terminal 
velocity on the effective temperature (Fig.|5J is more scattered. 
The variations of the mass-loss rate and the terminal velocity 
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M[M Q yr 1 ] (r eff =T effi0 ) v„ [km s~ 1 ] (T etr T eftfl ) 

Figure 5. Comparison of the mass-loss rates (left panel) and the terminal velocities (right panel) of studied stars calculated with original effective temperature 
and effective temperature higher by 1000 K. Solid line denotes one to one relation and dashed line denotes linear fit to the relation between two groups of 
wind parameters. 




M [M yr 1 ] (M=M ) v„ [km s~ 1 ] (M=M ) 

Figure 6. Comparison of the mass-loss rates (left panel) and the terminal velocities (right panel) of studied stars calculated with original mass and mass 1.05 
times higher. Solid line denotes one to one relation and dashed line denotes linear fit to the relation between two groups of wind parameters. 



with stellar mass are more stringent (Fig.|6j. With increasing stellar 
mass the mass loss rate decreases, on average 

M ~ M" 1 ' 60 , (7) 

with fair agreement with Eqs. J4j for the derived value a' — 
0.38. The terminal velocity is proportional to the escape velocity. 
From our calculations we derive the average relation 

Uoo ~M ' 52 , (8) 

which clearly reflects this proportionality. 

It is clear that some part of the discrepancy between observed 
and theoretical wind parameters may be attributed to the uncertain- 
ties in the determination of the stellar parameters (mass and the ef- 
fective temperature). Especially the downward revision of the stel- 
lar mass by the factor of 1.5x (e.g. due to the difference between 
spectroscopic and evolutionary masses) would cause downward re- 
vision of the terminal velocities by the factor of about 1.3 and the 
increase of the mass-loss rate by 2 x . On the other hand, since the 
distance to the SMC is known with a relatively high degree of preci- 
sion we conclude that probably also another source of discrepancy 
between observed and predicted wind terminal velocities is present. 
Moreover, if higher terminal velocities scatter is given purely by 



the uncertainties of determination of stellar mass and effective tem- 
perature, then similar effect should be also present in KK1, where 
slightly better agreement between theoretical and observed termi- 
nal velocities was derived. 

4.2 Abundances 

To study the variations of wind parameters with metallicity we 
recalculated wind models with the same stellar parameters as in 
Tab. [2] however with abundance of heavier elements 1.5 times 
higher. Comparison of the mass-loss rates and the terminal veloci- 
ties calculated with different metallicities is given in Fig.|JJ 

For higher metallicity the radiative force is higher and conse- 
quently also the mass loss rate is higher. We have found that for 
studied stars the relation 

M ~ Z 0S7 (9) 

holds. This is in a good agreement with VKL (see Eq. |3J)- 

The situation with the variations of the terminal velocity with 
the metallicity is more complicated. In Sect. l3.2l we concluded that 
the ratio of the wind terminal velocity to the stellar escape veloc- 
ity Voo/v^sc may be slightly higher for Galactic stars than for SMC 
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Figure 7. Comparison of the mass-loss rates (left panel) and the terminal velocities (right panel) of studied stars calculated with original metallicity and 
metallicity 1.5 times higher. Solid line denotes one to one relation and dashed line denotes linear fit to the relation between two groups of wind parameters. 
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Figure 8. The radial variation of the contribution of the radiative force /™ d 
acting on element h to the total radiative force / md for star AzV 238 (plotted 
as to™ d = /™ d // md ). Iron is the most important element for the radiative 
driving close to the star, while oxygen dominates in the outer regions. 



stars. With varying metallicity the wind density changes, conse- 
quently wind ionization and excitation state varies, and hence also 
terminal velocities of studied stars vary. However, these variations 
are not monotonic, as can be seen from Fig.0 For some stars the 
terminal velocity does not significantly change with metallicity, for 
some of them increases, for some of them decreases. On average, 
the terminal velocity slightly increases with increasing metallicity 
for studied SMC stars, 

Voo ~ Z 0M . (10) 
we can infer more steep pro- 



From calculations of Kudritzki 1200'. 
portionality, roughly v x 



7O.12 



4.3 Sensitivity of Vao on the conditions in the outer wind 

Although on average the value of depends on metallicity only 
slightly, the high scatter between terminal velocities derived using 
different metallicities in Fig.0may seem surprising. Similar effect 
was however, in a different context, reported also by other authors 
(e.g. Puis et al.l2000). The O star winds in the outer regions are ac- 
celerated mainly by few dozen lines of lighter elements (like C, N, 



O, see Pauldrach [l987l Vink et al. ll999l and also Figs.[S]and|9}. In 
such a case, the line acceleration is very sensitive to the detailed 
wind structure (i.e. temperature, density and chemical composi- 
tion). Partly due to this, large variations of Uoo/itsc (both observed 
and theoretical, see Lamers et al. 1 19951 and Pauldrach at al. ll990l) 
occur for individual stars. This is also likely cause for the high scat- 
ter between terminal velocities derived using different metallicities 
in Fig. This scatter also occurs when variations of terminal ve- 
locity with temperature (see Fig.|5J are studied. To a lesser extent 
this is also true for the variations of stellar mass, see Fig.|6| 

The metallicity of studied stars is not known with a high 
precision. For many stars we used just scaled Galactic chemical 
composition, although significant deviations from this scaling exist 
(e.g. B03; they are also probably manifested by a different number 
ratio of WN to WC stars in the Clouds, see Mikulasek 1969). More- 
over, we do not have any information on the abundance of many el- 
ements which are important for the radiative driving at all (e.g. Ne, 
Ar). Even worse, there are large differences in the chemical com- 
position of individual stars (e.g. B03). To conclude, some part of 
higher scatter between observed and predicted terminal velocities 
can be attributed to poorly known chemical composition of most of 
the studied stars. Uncertainties of the determination of other stellar 
parameters, i.e. the stellar mass and the effective temperature may 
(besides the approximations involved in our code) also contribute 
to this scatter. 



5 MULTICOMPONENT EFFECTS 

Stellar winds of hot stars are accelerated mainly by the absorption 
of radiation in the resonance lines of heavier elements. The radia- 
tive acceleration acting on individual elements is however differ- 
ent, and consequently individual elements have different velocities. 
Consequently, stellar winds of hot stars have a multicomponent na- 
ture (e.g. Springmann & Pauldrach 1 19921 KKII). For many stars 
the velocity differences are small, thus the wind multico mponent 
nature can be neglected in this case. Krticka et al. showed 
that in the low-metallicity environment the velocity differences be- 
tween wind components are larger, thus multicomponent effects are 
more important. However, as was noted e.g. by Mr04, even for rela- 
tively low-metallicity environment of stellar wind of SMC stars the 
multicomponent effects can be neglected. 
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Figure 9. Contribution of lines with different strengths to the radiative force at the critical point (r 
(r fa 18.5i?„, right panel) for star AzV 238. The total relative contribution to the radiative force /j 
contribution to the radiative force /k^£// is plotted using boxes. Solid line denotes number of spectral lines whose relative contribution to the radiative force 
lies in a given interval. Radiative acceleration close to the star is mainly due to hundreds of lines whose f^t/ f lies in the interval (10 -3 — 10~ 2 ), whereas 
the acceleration in the outer regions is given just by few dozens lines with 10~ 2 < /, rad // rad < 10 — 1 . 



On the other hand, these considerations are based on the as- 
sumption that the heavier elements can be described by one com- 
ponent only. By other words, usually it is assumed that atomic mass 
of heavier ions is the same, the radiative acceleration on heavier el- 
ements is the same and consequently these elements have the same 
mean velocity. However, in reality individual elements have dif- 
ferent mass, the radiative acceleration acting on them is also differ- 
ent and consequently these elements have different mean velocities. 
Hence, multicomponent effects may occur for higher densities and 
metallicities than previously assumed. 



5.1 Velocity differences in the stellar wind 

To estimate the importance of multicomponent effects in the stellar 
wind of SMC stars we calculate approximate velocity differences 
between individual heavier elements h (accelerated by the line ab- 
sorption) and passive wind component p (hydrogen and helium) 
using the momentum equation of individual wind elements. In the 
case of stationary spherically symmetric stellar wind this equation 
reads (cf. Burgers l 19691 KKII) 



dv r 



dr 



rad 
9a 



9 T [ a aPa) H E-\ 

p a dr m a 



Vrb — "ra 
\Vrb - Vra\ ' 



(ID 



where v ra and p a are the velocity and the density of component a 
(a = p or a — h), a a is the isothermal sound speed, E is the charge 
separation electric field, g is the gravitational acceleration, g a ad is 
the radiative acceleration and q a and m a are charge and mass of 
particle a. Frictional parameter has the following form: 



K ab = n a n b —— — In A, 

klab 



(12) 



where n a and rib are number densities of individual components 
and mean temperature of both components is 



Tab = 



m a T b + m b T a 
m a + rrib 



(13) 



(T a and T b are temperatures of wind components a and b) and 
In A is the Coulomb logarithm. The argument of the Chandrasekhar 
function G(x ab ) is 



\Vrb 



Ctab 



where 



2 2k (m a T b + m b T a ) 
a a b = - 



m a m b 



(14) 



(15) 



For low velocity differences (x ab < 1) the flow is well coupled. 
However, for higher velocity differences (x ab > 1) the Chan- 
drasekhar function is decreasing and this behaviour may enable dy- 
namical decoupling of wind components (see Springmann & Paul- 
drach ll99l KKII). 

In the momentum equation of heavier ions Eq. {TTJ the left- 
hand side term and the pressure term can be neglected. Also the 
gravitational acceleration and the electric polarisation field term 
can be neglected. Finally, due to high number density of passive 
component (hydrogen and helium) compared to the number den- 
sity of heavier elements the only important frictional term is that 
between passive and other components. Consequently, the approx- 
imative momentum equation of heavier element h is 



9h ad = — K hp G(x hp ) 

Ph 



(16) 



This equation states that whole momentum obtained by individual 
heavier elements due to the line-absorption is transferred by friction 
to the passive component. Using the Taylor expansion of the Chan- 
drasekhar function (appropriate for x^ v < 1) G(xhp) ~ §^7* the 
relative velocity difference between passive component and a given 
heavier ion h is 



Xhn — 



Vrh\ 



a In 



rad Ttlh 
9h 



3kT 



(17) 



where we have assumed T„ 



n p 80rg 2 (?^lnA' 
Th ~ T. For our discussion the rel 



evant quantity is not the radiative acceleration g}j ad itself, however 
the radiative force /£ ad = phg*h d (P er un it volume). Thus, 

1 3kT 



Xhp 



rrad 

Jh 



nhn p %^/nq\qp In A 



(18) 
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The work done by the radiative acceleration is used to lift the wind 
material from the stellar gravitational well. Thus, for stars with 
the same mass-loss rates (and the same velocity field) the radia- 
tive force (at corresponding radii) shall be similar. These stars with 
lower metallicities have lower rih and consequently higher Xhp- 
In reality, stars with lower metallicities have also lower mass-loss 
rates, lower n p and consequently even h igher Xh P - Due to these 
two effects (see also Krticka et al. 120031) stars with lower metal- 
licities have higher velocity differences between wind components. 
The crucial point is, however, that heavier elements which are not 
abundant in the stellar wind (i.e. iih <C n p ) and which significantly 
contribute to the radiative acceleration (i.e. their f™ d is large) may 
have large relative velocity differences Xhp, in many cases x/ lp m 1. 
In such a case, instability connected with the decoupling of consid- 
ered element may occur. 

Let us first roughly estimate in which situation this may occur. 
Let be tojj the relative contribution of a given element h to the 
total radiative force / rad , i.e. 



/rad 
h 



rad /-rad 
% J ■ 



(19) 



Neglecting the gravity, the pressure term and the electric polari- 
sation field the total radiative force can be approximated from the 
momentum equation Eq. 1 1 1 i of passive component as 



rrad _ "tV 



(20) 



where v r is the mean wind velocity, we may assume v r « v rp . The 
velocity gradient can be estimated as 

dv r l4 
^"IV (21) 

for rih we can write from the approximate continuity equation 



n h 



Z h 



M 



(22) 



where Zh is the density of a given element relative to the bulk den- 
sity p in the stellar atmosphere (ph = Zhp)- Consequently using 
Eqs. (HH22J we derive 



■l h n 



rad 3 D rn p m h 



Z h M In A 



(23) 



In the case when the heavier ions are approximated by one compo- 
nent only, we have w™ d = 1 and we arrive at Eq. (6) of Krticka 
et al. (2003). However, in many cases the abundance of individ- 
ual element Zh is much lower than 2) ; (the ratio of total density 
of heavier "absorbing" ions and passive component in the atmo- 
sphere) and thus, predicted velocity differences are larger. Eq. <23> 
can be also rewritten in a more convenient form as 

^. n ni c...rad..3 r> ^h 

<Khp 



0.015w T h d vlR 12 



Z h M. 



(24) 



where we have assumed m p = ran, q P = ez p , qh = ezh and 
mh = -AfcTOH, where e is the elementary charge and mu is the pro- 
ton mass, and scaled quantities are = M/(10~ 1 M0yr~ 1 ), 
1-8 = ^/(lO^ms- 1 ), R 12 = i?»/(10 12 cm), and T 4 = 
T'/(10 4 K). From Eq. 1241 it follows that for heavier elements 
(Ah ~ 10) with low abundance (Zh ~ 10 -5 ) which signifi- 
cantly contribute to the radiative force (i.e. Wh ~ 0.1) the de- 
coupling (Xh p > 1) can occur for relatively large mass-loss rates 
(of order 1O" 8 M yr _1 ). The mass-loss rate for which the decou- 
pling occurs in the winds of SMC stars is approximately five times 
higher than the mass-loss rate for which the decoupling occurs in 
the winds of Galactic stars, since Zh, smc ~ 0-2Zh,Q, 



Eq. 1241 provides only a very approximate expression for the 
velocity difference, mainly due to a very simplified velocity gra- 
dient assumed in Eq. <2U . For a more reliable estimate it is much 
better to use Eq. GHJ. 

Note that this picture of decoupling is somehow different from 
that presented by other studies, i.e. Springmann & Pauldrach { 19921) 
and KKII and discussed by Mr04. These studies assumed that all 
heavier ions which are accelerated by line-transitions decouple 
from the mean (passive) flow. However now we discuss a more re- 
alistic description that individual ions decouple from the mean flow 
separately. 



5.2 Velocity differences in the winds of studied SMC stars 

We used our NLTE wind code to calculate approximative velocity 
differences in studied SMC stellar winds. For this purpose we ap- 
plied Eq. 1 1 8i . where both the contribution of a given element to 
the radiative force Wh and charges of wind components qh and q p 
are calculated using our NLTE wind code. 

Calculated values of Xh p for individual stars and individual 
elements are given in Fig. 1101 Generally, the relative velocity dif- 
ferences are smallest close to the star, where the stellar wind is 
relatively dense. As the stellar wind accelerates, the wind density is 
lower and velocity differences are higher. At some point the veloc- 
ity differences have its maximum and for larger radii decrease out- 
wards due to decreasing velocity gradient. This general behaviour 
of relative velocity differences was described elsewhere (KKII). In 
our case the behaviour is more complicated mainly due to the pro- 
cesses of ionization and recombination. 

For most of the stars and for most of the elements the veloc- 
ity differences are rather low, Xh P <C 1, thus there occurs no de- 
coupling in this case. However, in some cases, especially in those 
where the wind density is low, the relative velocity difference be- 
tween argon (or sulphur) and passive component becomes higher 
and close to one. This is caused by the fact that argon and sulphur 
have very low metallicity (Zat, smc ~ 3 10~ 5 ) and significantly 
contribute to the radiative force (w™? « 0.1). The velocity differ- 
ences are so high, that they may cause large frictional heating. 



5.3 Multicomponent wind models 

To test the possibility of frictional heating, we have calculated four- 
component wind models, where the fourth component is either ar- 
gon, sulphur or carbon, depending on which element for a given 
star has a maximum velocity difference. The other three compo- 
nents are heavier ions, passive component (hydrogen and helium) 
and free electrons. 

First, we tested whether the simple equation Eq. 1 1 8i is able to 
reliably predict the velocity differences. We plot the result only for 
the case of star NGC 346 WB 1 (see Fig ll II . however the result that 
Eq. J I SI is able to very reliably predict the approximate velocity 
differences is similar also for other stars. 

For stars with highest velocity differences the frictional heat- 
ing occurs 1 . This is shown in Fies. 1121 1131 where we compare 



1 In the common case that occurs in the stellar atmosphere and wind it is 
usually not necessary to differentiate between kinetic temperatures of indi- 
vidual particles since these temperatures are nearly equal. However, in the 
present case the temperatures of individual wind components significantly 
differ. Since the collisional terms that enter into the statistical equilibrium 
equations are caused by the collisions with free electrons that move with 
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Figure 10. Relative velocity differences between individual heavier elements and passive component calculated using Eq. 1181 . Only values for elements with 
highest relative velocity differences are plotted. Velocity differences for other heavier elements are much lower. 
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Figure 12. Comparison of four-component model (with sulfur as a fourth component, dashed line) and three-component model (solid line) for star NGC 346 
MPG 487 (left) and NGC 346 MPG 113 (right). Top panel: Comparison of velocity of passive component (hydrogen+helium) in individual models. Note that 
the velocities of all individual wind components in four-component or three-component models are similar. Middle panel: Comparison of wind temperatures. 
Wind temperatures of individual wind components in three-component model are nearly equal. Wind temperature in four-component model is higher due to 
the frictional heating and sulfur temperature is even higher than the temperature of other components. Bottom panel: Relative velocity difference between 
passive component and sulphur in four-component model. The velocity difference in the case of NGC 346 MPG 1 1 3 is close to one in the outer wind regions, 
thus potentially enabling instability due to sulphur decoupling. 



four-component models (with wind components sulphur, remain- 
ing heavier ions, passive component and electrons) with three- 
component models (with wind components heavier ions, passive 
component and electrons). In the outer wind regions the relative 

much higher thermal speed than other particles, we inserted the electron 
temperature into all terms in the statistical equilibrium equations. This is in- 
deed only approximative in the case of relatively large velocity differences 
between wind components. 



velocity difference between passive component and sulphur x f s is 
close to 1. Due to this higher velocity difference the stellar wind is 
frictionally heated in the outer parts. Consequently, wind temper- 
ature is higher in four-component model than in three-component 
one. Higher wind temperature also slightly modifies the radiative 
force since wind velocity is slightly higher in four-component case. 

The frictional heating also slightly modifies the ionization 
equilibrium. The ionization fraction of higher ionization stages 
(such as Nv, Ov) of frictionally heated wind is higher than in 
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Figure 13. The same as Fig. ll2l for SMC-N81 stars #3 (left) and #1 1 (right). Note that for winds of both stars the friction between sulphur and passive wind 
component (hydrogen and helium) is important for wind structure. In both cases the sulphur decoupling may occur. In the case of SMC-N8 1 #3 sulphur may 
decouple for velocities lower than the escape speed (the surface escape speed for this star is v esc = 1180 km s -1 ) and in the case of SMC-N81 #1 1 sulphur 
may decouple for velocities slightly higher than the escape speed (the escape speed is also plotted on the graph of this star). 



the wind with negligible frictional heating, whereas the ionization 
fractions of lower ionization stages (e.g. C IV, Si IV) are lower in 
the frictionally heated wind. The differencies are by about a factor 
of 2. 

The velocity differences in the case of frictionally heated four- 
component wind models are larger than those estimated using ap- 
proximate Eq. 1 1 SI (cf. Fig. llOI . since velocity difference increases 
with increasing temperature. The velocity difference in the case of 
NGC 346 MPG 113, SMC-N81 #1 1 and #3 is so high that the wind 
instability connected with sulphur decoupling may occur (Owocki 
& Puls l2002l Krticka & Kubat|2002). For star #3 this occurs for 
velocities lower than the escape speed. 



We have shown that for stars that have very low observed 
mass-loss rates the velocity differences for some heavier elements 
are high. In such case either frictional heating becomes important or 
even some elements may decouple from the wind. Note that these 
multicomponent effects are important already for theoretically de- 
rived mass-loss rates, i.e. those that are much higher than the ob- 
served ones. Since the frictional heating increases the wind tem- 
perature only in the outer wind regions (for velocities higher than 
that corresponding to the critical point velocity), its influence on 
the mass-loss rate is negligible. 

Fall back of the wind material (which may cause the lowering 
of wind mass-loss rate) may occur only if the wind components de- 
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Figure 11. Comparison of relative velocity differences between passive 
component and argon calculated using approximate expression Eg. 1181 and 
using four-component model. Clearly, in the case of negligible frictional 
heating Eg. 1181 is able to correctly predict velocity differences. 



couple below the point where wind velocity is equal to the escape 
velocity. Only in this case the runaway instability connected with 
decoupling may lead to fall ba ck of wind material (Porter & Sk- 
ouza[l999, Votruba et al. 2006) and consequently lower the mass- 
loss rate. This occurs only for star SMC-N81 #3. For other stars we 
have not obtained the possibility of decoupling close to the stellar 
surface, for velocities lower than the escape velocity. This means 
that presented models are, with one exception, close to the star sta- 
ble for perturbations of multicomponent flow. Consequently, we are 
not able to explain very low observed mass-loss rates using multi- 
component models for all except one star. Only for star SMC-N81 
#3 the low observed mass-loss rate may be connected with multi- 
component effects. 



6 DISCUSSION AND CONCLUSIONS 

We have presented NLTE wind models of cooler O stars in SMC. 
These models enable prediction of basic wind parameters, i.e. the 
mass-loss rates and the terminal velocities. We have compared 
these predicted wind parameters with those derived from obser- 
vation. We have concluded that there is a relatively good agree- 
ment between observed and predicted terminal velocities with some 
scatter. This scatter can be partly attributed to the poorly known 
chemical composition of SMC stars. Moreover, the existence of 
unidentified binaries in our sample can also influence the com- 
parison. Last but not least, the simplification involved in our code 
(e.g. simplified treatment of radiative transfer equation, like neglect 
of line overlaps) can also influence the results. However, the situa- 
tion with mass-loss rates is different. For relatively high mass-loss 
rates (i.e. M > lO _7 M0yr _1 ) there is again a good agreement 
between predicted quantities and those derived from observations. 
However, for smaller mass-loss rates the agreement is rather poor, 
predicted mass-loss rates are significantly higher than those de- 
rived from observation. Note that this is not an effect of our models 
only, since there is a similar disagreement between observations 
and models of VKL (see B03 and Mr04). 

We have tested whether the mentioned systematic disagree- 
ment between theoretical and observed values of mass-loss rates 
can be caused by the decoupling of some metals from the mean 
flow. We have shown that for stars with relatively good agreement 



between observed and predicted mass-loss rates the velocity differ- 
ences are small (and thus there is no decoupling of wind compo- 
nents). On the other hand for stars for which the predicted mass- 
loss rates are systematically too high either the instability con- 
nected with the decoupling of wind components may occur or the 
frictional heating increases wind temperature. Note however that 
only for one star we have obtained the possibility of decoupling 
close to the stellar surface for velocities lower than the escape ve- 
locity. (Note that in order to obtain fall back of the wind material 
and consequently lower wind mass-loss rate it is necessary to the 
decoupling occur below the point where wind velocity is equal to 
the escape velocity.) With this exception, presented models are, to 
our knowledge, close to the star stable for perturbations of mul- 
ticomponent flow. This means that we are not able without any 
further assumptions to explain very low observed mass-loss rates 
of most studied stars with M < M@ yr" 1 (although in the 

outer wind regions of these stars the multicomponent effects are 
important). However, there are several possibilities how to (from 
the theoretical point of view) obtain decoupling of wind compo- 
nents for the velocities lower than the escape velocity also for other 
stars (like slightly different metallicity, overestimation of theoreti- 
cal mass-loss rates, neglect of some heavier elements or improved 
treatment of the multicomponent flow). Another way how to ex- 
plain low wind mass-loss rates may be connected with the thin- 
wind effect (Puis et al. ll998l Owocki & Puls ll999l) . 

The decoupling discussed here slig htly dif fers from that stud- 
ied so far (cf. Springmann & Pauldrach |l992j, KKII). In the pre- 
vious studies it was usually assumed that metals and passive com- 
ponent may decouple completely. We present here a more detailed 
and more realistic picture of the decoupling of wind components. 
We showed that some elements which have very low abundance, 
however which significantly contribute to the radiative force, may 
decouple (or cause the instability connected with the decoupling) 
individually from the mean flow, that is basically not decoupled. 
Moreover, the decoupling may occur for higher mass-loss rates than 
it was previously assumed. We present an approximate formula for 
the test of importance of discussed decoupling of individual heavier 
elements. 

For our comparison we used mass-loss rates derived from 
the observation assuming smooth winds. If the O star winds are 
clumped, as is indicated by observations of e.g. B03 or Martins et 
al. 120051) . then the mass-loss rates of O stars are likely lower than 
predicted by current hot star wind theory and we obtain a disagree- 
ment between observations and theory for all stars. On the other 
hand a better treatment of the opacity sources in the UV domain 
with account of the line overlaps and comoving-frame radiative 
transport may help to overcome this potential disagreement. 

We have studied the dependence of wind parameters on metal- 
licity. It is well-known that the mass-loss rate increases with in- 
creasing metallicity. We have derived that the mass-loss rate scales 
with metallicity as M ~ Z ' 67 . The terminal velocity of individ- 
ual stars also varies with metallicity mainly due to the sensitivity 
of the radiative force on detailed wind state in outer wind regions. 
However, on the average the terminal velocity varies with metallic- 
ity only slightly (for studied values of metallicity) as v x ~ Z 06 . 
As a consequence, the ratio of wind terminal velocity to the sur- 
face escape velocity Hoo/^csc ~ 2.3 is only slightly lower than for 
Galactic stars. 
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